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Abstract 

We study the properties of a one-dimensional (ID) granular gas consisting of N hard rods on 
a line of length L (with periodic boundary conditions). The particles collide inelastically and 
are fluidized by a heat bath at temperature Tb and viscosity 7. The analysis is supported by 
molecular dynamics simulations. The average properties of the system are first discussed, focusing 
on the relations between granular temperature T g = m(v 2 }, kinetic pressure and density p = N/L. 
Thereafter, we consider the fluctuations around the average behavior obtaining a slightly non- 
Gaussian behavior of the velocity distributions and a spatially correlated velocity field; the density 
field displays clustering: this is reflected in the structure factor which has a peak in the k ~ 
region suggesting an analogy between inelastic hard core interactions and an effective attractive 
potential. Finally, we study the transport properties, showing the typical sub-diffusive behavior 
of ID stochastically driven systems, i.e. {\x(t) — x(0)| 2 ) ~ Dt 1 ^ 2 where D for the inelastic fluid is 
larger than the elastic case. This is directly related to the peak of the structure factor at small 
wave- vectors. 

PACS numbers: 02.50.Ey, 05.20.Dd, 81.05.Rm 



I. INTRODUCTION 



Scientists and engineers have been studying granular materials for nearly two centuries 
for their relevance both in natural processes (landslides, dunes, Saturn rings) and in indus- 
try (handling of cereals and minerals, fabrication of pharmaceuticals etc.) . 1 ' 2 ! 3 ' The under- 
standing of the "granular state" still represents an open challenge and one of the most active 
research topics in non-equilibrium statistical mechanics and fluid dynamics. For instance, a 
way to attack the problem consists in fluidizing the grains by shaking them so that the sys- 
tem behaves as a non-ideal gas, a problem relatively easier to study. The difficulty, but also 
the beauty, of the dynamics of granular gases, meant as rarefied assemblies of macroscopic 
particles, stems from the inelastic nature of their collisions which leads to a variety of very 
peculiar phenomena. Several theoretical methods have been employed to deal with granular 
gases ranging from hydrodynamic equations, kinetic theories to molecular dynamics. Engi- 
neers often prefer the strategy of the continuum description because it gives a better grasp 
of real life phenomena, while natural scientists tend to opt for a microscopic approach, to 
better control each step of the modelization. The latter, as far as the interaction between 
particles is concerned, regards granular systems as peculiar fluids, and treats them through 
the same methods which have been successfully applied to ordinary fluids. 5 This allows not 
only, to employ concepts already developed by physicists and chemists, but also to stress 
analogies and substantial differences. 

The purpose of the present paper is to establish such a connection for a system of stochas- 
tically driven inelastic hard-rods constrained to move on a ring. The elastic version of this 
system has a long traditio n 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 and is particularly suitable to test approximations 
and theories since many of its equilibrium properties can be derived in a closed analytical 
form. Even though, the one dimensional geometry introduces some peculiarities not shared 
by real fluids, we shall show that the model provides many useful information and a very 
rich phenomenology which closely recalls the behavior of microscopic particles confined in 
tubules or cylindrical pores with little interconnection. A second reason to investigate such 
a model is to show how the inelasticity of interactions influences not only the average global 
properties of a system, but also its microscopic local structure. 

A basic requirement to a theoretical description of a granular gas is to provide an equation 
of state linking the relevant control parameters and possibly to relate it to the microscopic 
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structure of the system. This connection is well known for classical fluids, where thermo- 
dynamic and transport properties are linked to the microscopic level via the correlation 
function formalism. 

One dimensional models have been employed by several authors as simple models of 
granular gases . 13 1 14 1 15 1 16 i 17 i 18 1 19 1 20 i 21 1 22 1 23 i 24 The differences between the various models stem 
chiefly from the choice of the thermalizing device. In fact, granular gases would come to 
rest unless supplying energy compensating the losses due to the inelastic collisions. We 
call, by analogy, "heat bath" the external driving mechanism maintaining the system in a 
statistically steady state. 

For the history, the first ID models, which were proposed, had no periodic boundary 
conditions and the energy was injected by a vibrating wall (stochastic or not). This kind of 
external driving however, was not able to keep the system homogeneous, because only the 
first and last particle had a direct interaction with the wall.— As an alternative, a uniform 
heating mechanism, namely a Gaussian white noise acting on each particle, was introduced.— 
Later Puglisi et al.— added a second ingredient, consisting of a friction term that prevents 
the kinetic energy from diverging. With such a modification the system reaches a steady 
regime and time averages can be safely computed. In the present paper we shall focus on 
this last model characterizing its steady state properties. 

The layout is the following: in section II we introduce the model, in section III we obtain 
numerically and by approximate analytical arguments equations for the average kinetic 
energy and pressure. Section IV is devoted to fluctuations of the system observables around 
their average values. In section V we study the diffusion properties of the system. Finally, 
in section VI we present the conclusions. 

II. THE MODEL 

Inelastic hard sphere models are perhaps the simplest models able to capture the two 
salient features of granular fluids, namely the hard core repulsion between grains and the 
dissipation of kinetic energy due to the inelastic collisions. Since many of the equilibrium 
properties of the ID elastic hard rods are known in closed analytical form, such a system 
represents an excellent reference model even for the inelastic case. Let us consider N identi- 
cal impenetrable rods, of coordinates Xi(t), mass m and size a, constrained to move along a 
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line of length L. Periodic boundary conditions are assumed. The hard-core character of the 
repulsive forces among particles reduces the interactions to single binary, instantaneous col- 
lision events occurring whenever two consecutive rods reach a distance di(t) = Xi(t) —Xi-i(t) 
equal to their length a. When two inelastic hard-rods collide, their post-collisional velocities 
(primed symbols) are related to pre-collisional velocities (unprimed symbols) through the 
rule: 

1 + r 

v'i=Vi — (Vi-Vj) (1) 

where r indicates the coefficient of restitution. The interaction of each particle with the 
heat-bath is represented by the combination of a viscous force proportional to the velocity 
and a stochastic force. Then each particle follows the so called Kramers dynamics 

m ~dt = ~ m ^ Vi + ^ t ) ( 3 ) 

where 7 is the viscous friction coefficient, is a Gaussian white noise with zero average 
and correlation 

m)^)) = ^rnT b 6 t3 6(t-s) , (4) 

T;, is the "heat-bath temperature" and (■) indicates the average over a statistical ensemble 
of noise realizations. 

We have developed a numerical simulation code for hard rods interacting through mo- 
mentum conserving but energy dissipating collisions. In our simulations the motion between 
two consecutive collisions is governed by the dynamics (|2I3|) . Thus, we determine the instant 
when the first collision among the N particles occurs and change their velocities and posi- 
tions according to the equations of motion. The effect of the collision is taken into account 
by updating the velocities after each collision according to the rule (JTJ. 

We tested our code on the elastic case (r = 1) and checked that our simulations faithfully 
reproduced the well known properties of the equilibrium hard rod system. 

III. AVERAGE PROPERTIES 

We begin by considering the steady state properties of the model. The aim is to derive 
relations connecting the microscopic parameters to the "thermodynamic" observables such 
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as temperature and pressure and eventually to obtain an "equation of state" relating these 
two quantities. In order to achieve this goal we assume that the system is homogeneous, so 
that its density p is constant. 



A. Kinetic temperature 

Collisions and the Kramers' dynamics entail that the time derivative of the average kinetic 
energy per particle is 

j t 1 -m(v(tf)= 1 (T b -T g )-w(t) } (5) 

where T g = m^2f =1 vf /N is the granular temperature and w(t) is the average power dis- 
sipated by collisions, given by w = ^—^-m{5v 2 ) /r c , where 5v is the difference between the 
pre-collisional velocities of the colliding pair. The average collision time r c is estimated by 
assuming a mean free path A = (L — aN)/N, where L — oN is the free volume. We obtain, 
in terms of the system density p = N/ L, 

A 1 - op fm 

V P V Tg 

Thus the average power dissipated per grain reads: 



W =—(5v) T -—. (7) 

In order to estimate T g we assume that m(5v 2 ) ~ AT g since the velocities of the colliding 
pairs are strongly correlated. Thus imposing the solution of Eq. (JHJ) to be stationary we 
obtain for T g the following expression: 

Tg = !=■ (8) 

\ _i_ 1-r 2 p T g 
27 1— per y m 

In figure ^ we compare formula (JHJ with the results of numerical simulations at various 
densities. In spite of the simplicity of the argument used to derive Eq. (jHJ), the agreement 
between T g extracted by simulations and its theoretical estimate is fairly good. 



B. Kinetic pressure 

In a granular system the total pressure, P, can be obtained via its mechanical or kinetic 
definition, i.e. as the impulse transferred across a surface in the unit of time . 25 ' 26 The pressure 
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FIG. 1: Comparison between the numerical results for the granular temperature versus density and 
the corresponding theoretical expression Eq. (|5|>. The simulation data refer to a = 0.2, 7 = 0.2, 
T;, = 1.0, r = 0.6 and r = 0.8. We kept the system size fixed to L = 40 but varied the number of 
particles N to change the density. 

contains both the ideal gas and the collisional contribution (P id and P exc respectively) 

P = P id + P exc = P T g {p) + -f - V 5 Pk (9) 

where the second equality stems from the virial theorem.^ 7 - Here t Q b is the observation time, 
the sum runs over the M c collisions and 5p k = mSvk represents the impulse variation due to 
the k-th collision. 

An approximate formula for P exc can be derived as follows. The average collision fre- 
quency per particle can be estimated as r" 1 = (M c /t b)/N. By replacing in Eq. @ £ & with 
(M c /N)t c and using r c given by Eq. (jUJ), we obtain, for the excess part of the pressure, the 
expression 

Pe XC = ^m{5v c ) = -^-T g {p). (10) 
t c L 1 — ap 

Collecting pieces together we arrive at 
P{p)=T g {p) 



A a P 



TM^~ (ID 



1 — ap 

which reproduces the well known Tonks formula 7 in the case of elastic particles and consti- 
tutes the sought equation of state for the inelastic system. Let us recall that in the elastic 
case, Eq. (jllj) can be written in the virial form 

P(p)=T g (p)p[l+ P ag(a)) (12) 



showing the connection between the macroscopic and the microscopic level, since g(a) = 
1/(1 — pa) is the equilibrium pair correlation at contact. 

We see, from figure (J2J), that the presence of the pref actor T g (p), which is decreasing 
function of the density, makes P(p) to increase more slowly than the corresponding pressure 
of the elastic system in the same physical conditions (i.e. same density and contact with the 
same heat bath). 

Equations (jSJ) and (fTUj) for temperature and pressure coincide, in the limit 7 — > 0, a — > 0, 
7T{, = Q = const, with those derived by Williams and MacKintosh^ .. 

Following the standard approach to fluids, we define, even for the inelastic system, the 
response of the density to a uniform change of the pressure for a fixed value of the heat bath 
temperature: 

XT = -we (13) 



which is plotted in figure E3 We observe that the response of the inelastic system to a 
compression is much larger than the corresponding elastic system at the same density, due 
to the tendency to cluster. 



2.5 



P/P 



2.0 

id 
1.5 



L.O 



0.5 



Ho 



/ 

/ 

/ 

/ 

— / 


/I 1 1 1 1 


1 ' ' 


/ 


— r= 1.0 






• r = 0.8 






* r = 0.6 

— theory 












[ , T T 1 


1 



0.2 



0.4 



pa 



0.6 



FIG. 2: Comparison between the kinetic pressure obtained from the simulations and the prediction 
of Eq. ^2 a t different densities for inelastic hard rods with coefficients of restitution r = 0.6 and 
r = 0.8. The remaining parameters are the same as in figure ^ Dashed line refers to the pressure 
of the corresponding elastic system. 
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FIG. 3: Compressibility, computed from Eq. ((T3]) . plotted versus density, for inelastic hard rods 
with coefficient r = 0.6 (the remaining parameters are the same as in figure Symbols refers 
to simulations while solid lines are the theoretical predictions obtained via the expression for the 
pressure (jllj) . The curve for the elastic system r = 1 (dashed line) is also reported for sake of 
comparison. In the inset shows the ratio between the inelastic and elastic compressibility. 

IV. FLUCTUATIONS 

So far we have considered only the global average properties of the granular gas. It is 
well known, on the other hand, that these system may exhibit strong spontaneous deviations 
from their uniform state. In this section we shall study fluctuations of the main observables 
in order to understand the qualitative effect of inelasticity on such a peculiar fluid. 

A. Velocity distributions 

One of the signatures of the inelasticity of the collisions is represented by the shape of the 
velocity distribution function (VDF), P(v). Non-Gaussian VDF's, displaying low velocity 
and high velocity overpopulated regions, have been measured experimentall y 28 i 29 i 30 i 31 i 32 i 33 
and in numerical simulations . 21 > 34 In figure|3]we show two VDF corresponding to two different 
values of r. 

Theoretical, numerical and experimental studies have shown that the VDF for inelastic 
(r < 1) gases usually displays overpopulated tails. The literature seems to indicate the 
lack of a universal VDF: in d > 1 the solution of the homogeneous Boltzmann equation 
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FIG. 4: Velocity distributions P(i>) for three values of the coefficient of restitution r = 1.0 (circles), 
0.8 (squares), 0.6 (triangles). The remaining parameters are N = 1000, L = 1000, a = 0.2, 7], = 1.0, 
7 = 0.2. Dashed lines indicates the corresponding Gaussian fit. 

with inelastic collisions (with a stochastic driving similar to ours but without viscosity) has 
overpopulated tails of the kind ~ exp(— v 3 / 2 ). 35 



B. Energy fluctuations 

Interestingly, the energy fluctuations of our system, E = m vf/2, display a scaling with 
respect to the number of particles. We are interested in the quantity ((E 2 ) — (E) 2 )/ (E) 2 as 
a function of iV (at fixed density p). Defining (v n ) = j v n P(v)dv, we have 

(E 2 ) - (E) 2 = ^J2 ( W> - <«?><«?>) ( 14 ) 
Under the hypothesis that the variables Vi are independently distributed, we get 

(E 2 )-(E) 2 = ^[N(v*)-N(v 2 ) 2 ] (15) 



Since for Gaussian variables (v A ) = 3(v 2 ) 2 , we find 

(E 2 ) - (E) 2 2 



(16) 



{E) 2 N 

which is a well known formula for equilibrium systems^. This scaling is fairly well verified 
in figure El 

In the case of a granular fluid, P{vj) is no more Gaussian, exhibiting fatter tails, so 
one observes (t> 4 ) > 3(t> 2 ) 2 (for example in the case r = 0.6 we have (v 4 )/(v 2 ) ~ 3.3). 
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FIG. 5: Energy fluctuation as a function of the number of particles, N for r = 1.0,0.8,0.6. The 
elastic case agrees with the theoretical prediction 2/N, whereas the inelastic case gives a value of 
the relative fluctuation slightly larger. The remaining parameters are the same as in figure |1] 

This leads to the conclusion that the scaling ~ 1/N of formula (|16|) still holds, but with a 
coefficient larger than 2. Simulation runs for r = 0.6 confirm this prediction (Fig. EJ). The 
renormalization of the multiplicative constant occurring in the inelastic system could be 
interpreted also as an "effective reduction" of the number of degrees of freedom. Indeed, the 
inelastic system has the tendency to cluster, as it will be shown, and therefore the effective 
number of independent "particles" appears smaller. Another appealing interpretation is 
that the inelastic systems possesses an effective "specific heat" larger than that of elastic 
systems. 

C. Velocity correlations 

A universal signature of the inelasticity is the presence of correlations between the ve- 
locities of the particles. We measured the structure function of the velocities Vi, S v (k) = 
(v(k)v(— k)), where v(k) is the Fourier transform of t>j. In figure El we show three S v (k) 
corresponding to the elastic (r = 1) and inelastic system with r = 0.8 and r = 0.6. As men- 
tioned above, the elastic systems is characterized by uncorrelated velocities and this reflects 
on a constant structure function. A certain degree of correlation is instead evident in the 
inelastic system. In fact, the inelasticity reduces by a factor r the relative velocity of two 
colliding particles and this leads to an increasing correlation among velocities. However, the 
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FIG. 6: Structure function of the velocity field for elastic and inelastic systems. The control 
parameters are the same as in figure 0] 

noise induced by the bath competes with these correlations, making the structure function 
not very steep. More specifically, S v (k) can be fitted, in the middle range of k values, by 
an inverse power ~ k~ 5 , while at high k values it reaches a constant plateau. This is the 
fingerprint of a persistent internal noise (velocity fluctuations are not completely frozen by 
inelastic collisions). 37 

D. Distribution of interparticle spacing and of collision times 

The probability distribution, P(Sx), of distances between nearest neighbor particles 5x = 
Xi — Xi-i, shown in figure provides information about the spatial arrangement of the 
system. In the elastic case one easily finds 



for Sx > a and for 5x < a with A = (1 — pcr)/p. The presence of inelasticity modifies such 
a simple exponential law in the way shown in figure In this case, the probability of finding 
two particles at small separation increases together with that of finding large voids. Such 
a picture is consistent with the idea of the clustering phenomenon:— two particles, after 
the inelastic collision, have a smaller relative velocity and therefore reach smaller distances, 
eventually producing dense clusters and leaving larger empty regions with respect to the 
elastic case. 




(17) 
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FIG. 7: Distributions of distances between nearest neighbor particles 5x = Xj — for elastic and 
inelastic systems, with the same parameters used before. The solid line indicates the exponential 
expected in the elastic case (see text). The state parameters are the same as in Fig. @] 

On the contrary, the probability distribution of collision times, shown in figure |Hl appears 
to always follow the theoretical (elastic) form P(t) = l/r c exp(— t/r c ). Apart from a trivial 
rescaling due to the change of the thermal velocity with r, it seems not to depend appreciably 
on the coefficient of restitution. 
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FIG. 8: Distributions of collision times. The dashed line indicates the exponential law expected 
for the elastic system (see text). Control parameters as in figure [Tj 

Such a finding is in contrast with the situation observed in 2D vibrated granular sys- 
tem a 30 i 39 . A possible explanation for this discrepancy is the following: in the inelastic ID 
system, there is correlation between the relative velocities and the free-paths (or free-times), 
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otherwise the distribution P(t) and P(Sx) would have had the same shape due to the trivial 
relation x = vt, v being the the average velocity of the rods. In particular, the fact that the 
peak of P(Sx) in 5x = does not yield a corresponding peak in the t = region of P(t) 
suggests that the shorter the distance between particles the smaller their relative velocity. 



E. Density fluctuations 

We turn, now, to the study of the structural properties of the inelastic hard-rod gas, by 
considering the pair correlation function or the structure factor. As mentioned in section II, 
the virial equation relates the pressure of an elastic hard-rod system to its microscopic 
structure. In the presence of inelasticity, however, we expect that the tendency to cluster is 
mirrored by a change in the structural properties of the fluid. Therefore we considered the 
behavior of the static (truly speaking steady state) structure factor 

S(k) = 1 + p J dx[g{x) - l}e ikx (18) 

for different values of pa and inelasticity. 

The spatial structure of the system is determined, as in ordinary fluids, by the strong 
repulsive forces. Their role is seen in the oscillating structure of g(x). The inelastic nature of 
the collisions provides a correction to g(x), which can be better appreciated by studying the 
small wavelength behavior of S(k) which develops a peak at small k recalling the Ornstein- 
Zernike behavior 

sw * s-W+o*' ' (19) 

The coefficient C2 is negative for hard rods, whereas it is positive for the inelastic system. 
For hard rods, S(k) is known^ and reads: 

S ( k ) = —77—-— ^77—77- (20) 

with b = -r^—. 

l—pa 

Figure El shows the typical behaviors of S(k) for elastic and inelastic systems. The numer- 
ically computed structure factor of the elastic system agrees rather well with equation (|2()j) . 
The inelastic system, instead, displays a peak in the small k region reflecting the tendency 
of the fluid to cluster. The peak increases with inelasticity, demonstrating that the energy 
dissipation in collisions is responsible for these long range correlations. Incidentally, we 
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FIG. 9: Typical behavior of the structure function of S(k) for different values of the coefficients 
r = 1.0,0.8,0.6. Notice the growth of the peak at small k when r decreases. 

comment that such a behavior of S(k) could be attributed to the presence of a long range 
attractive effective potential between the rods, as a result of dynamical correlations. 40 

Mackintosh and Williams 20 found that, in the case of randomly kicked rods in the absence 
of viscosity, the pair correlation function decays as an inverse power law, g(x) oc x~ v with 
?7 — > for r — > 1 and rj — > 1/2 for r — > 0. Correspondingly, one expects S(k) to diverge as 
k~ l l 2 as k — > 0, for very inelastic systems. In other words, the inelasticity leads to long range 
spatial correlations which are revealed by the peak at small k of S(k). We remark that, in 
spite of the apparent similarity between the equations of state for elastic and inelastic system, 
their structural properties are radically different. Such a phenomenon is the result of the 
coupling of the long- wavelength modes of the velocity field with the stochastic non-conserved 
driving force. In fact, due to the inelastic collisions, the velocities of the particles tend to 
align, thus reducing the energy dissipation. On the other hand, these modes adsorb energy 
from the heat bath and grow in amplitude, and only the presence of friction prevents these 
excitations from becoming unstable. The density field, which is coupled to the velocity field 
by the continuity equation, also develops long range correlations, and the structure factor 
displays a peak at small wave-vectors. 
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V. TRANSPORT PROPERTIES 



One dimensional hard-core fluids exhibit an interesting connection between the micro- 
scopic structural properties and diffusive ones. In the present section, we present numerical 
results for the collective diffusion and for the diffusion of a tagged particle and show how 
these are connected to the structure. 

A. Collective Diffusion and Self-Diffusion 

Let us turn to analyze the perhaps simplest transport property of the hard rods system, 
namely the self-diffusion, i.e. the dynamics of a grain in the presence of iV — 1 partners. 
The problem is highly non-trivial since the single grain degrees of freedom are coupled to 
those of the remaining grains. Such a single-filing diffusion is also relevant in the study of 
transport of particles in narrow pores.— 

The diffusing particles can never pass each other. The excluded volume effect represents 
a severe hindrance for the particles to diffuse. In fact, a given particle in order to move must 
wait for a collective rearrangement of the entire system. Only when the cage of a particle 
expands, the tagged particle is free to diffuse further. This is a peculiar form of the so called 
cage effect which is enhanced by the one-dimensional geometry. In addition, the cage effect 
produces a negative region and a slow tail in the velocity autocorrelation function. 

As an appropriate measure of the self-diffusion, we consider the average square displace- 
ment of each particle from its position at a certain time, that we assume to be t = without 
loss of generality 

m = Ej^frW -*«»]■) (21) 

At an early stage, the self- diffusion is expected to display ballistic behavior, R(t) ~ M 2 £ 2 , 
with \v | 2 = Tg/m, before any perturbation (heat bath and collisions) change the free motion 
of particles, i.e. when t <C I/7 and t <C r c . 

A system of non-interacting (i.e. non-colliding) particles subjected to Kramers' dynam- 
ics (I2I3|) displays, after the ballistic transient, normal self- diffusion of the form R(t) ~ 2D t 
with Do = Tg/j. This is well verified in figure ITU1 (circles) . 

Lebowitz and Percus^ 2 - studied the tagged particle diffusion problem for systems governed 
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by non-dissipative dynamics without heat bath and found a diffusive behavior described by 

R(t) = 2D cM t = X(\v\)t (22) 
with A = (1 — per) I p. Since = y^T/2-Km one obtains 

Dcou = — — —\U (23) 

p V 2nm v ; 

On the other hand, almost in the same years, Harris^ studied the behavior of R(t) in the 

case of N identical Brownian particles with hard-core interactions , i.e. obeying a single- filing 

condition, and obtained a sub-diffusive behavior increasing as 

R{t) = 2X^fj' 2 (24) 

where D is the single (noninteracting) particle diffusion coefficient. 

In figure EH we study the self-diffusion R(t) for elastic and inelastic systems in the 
presence of heat bath and viscosity, obtaining two different regimes separated by a typical 
time t c . In the first transient regime t < t c we we observe the ballistic motion. In the 
second stage, instead, we expect the sub-diffusive behavior, R(t) ~ t 1//2 , predicted by Harris 
and other authors . 4 ? 144 ' 4 ^ . The inelastic system displays the same sub-diffusive behavior, but 
with a multiplicative constant larger than 1, i.e. at equal times the granular (inelastic) fluid 
has a larger absolute value of R(t). 

It is interesting to analyze the connection between this transport property and the com- 
pressibility of the system, as remarked by Kollmann.— 

B. Connection between self-diffusion and structure 

We follow Alexander and Pincus 4 ^ argument in order to derive a formula for the self- 
diffusion and show the connection with the compressibility of the system. Let us consider 
the two time correlator: 

(Pk(t)p-k(0)) = J2(e lkxdt) e- lkx ^) (25) 

ij 

we set Xi(t) = Xi + (xi(t) — X{) = Xi + Ui, where Xj are the nodes of the ID lattice Xj = ia, 
(a being the lattice spacing). Expanding around Xj 

(p fc (t)p_ fc (0)> ~ k 2 ^2e ik ^- x ^( Ui (0) Uj (t)) (26) 

ij 
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FIG. 10: Self-diffusion: behavior of R(t) for three different systems: one without collisions (free 
particles), one with elastic collisions and the third with inelastic collisions (r = 0.6). 



hence 



Nk 2 

We assume, now, that the density correlator varies as 



Ufc(t)«_fc(0)) = — 2 (27) 



(Pk(t)p-M) = {Pk(0)p-k(0))e- v ^ (28) 

where D = D /S(0) (this is demonstrated in Appendix A) is the collective diffusion coeffi- 
cient, with Dq = Tf,/j. Now, the mean square displacement per particle (}2T|) can be written 
as R(t) = ^2i([ui(t) — ui(0)] 2 )/N and, in Fourier components, reads 

R(t) = y £<u fc (0)u_ fc (0) - u fc (*)u_ fc (0)) (29) 



N 
k 



employing Eq. 123 we find 



j _ e -Dk 2 t 

R(t)=2j2(pk(0)P-k(0)) (30) 

k 

Approximating the sum with an integral (Ylk L/(2n) j^ a ,dk and recalling that 

NS(k) = (p fc (0)/j_ fc (0)> (31) 

we obtain 

„, . 2L f n/a dk n „A -e~ Dk2t 

m =Nj (32) 
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and therefore 

R(t) ~ -S(0)VM = - JMM (33) 
pn p V 7r 

Notice that such a formula in the case of hard- rods is identical to formula (J2l|) . 

We see that the tagged particle diffusion depends on the structure of the fluid. In the 
granular fluid the k — > part of the spectrum is enhanced and thus we expect a stronger 
tagged particle diffusion. This is what we observe. Physically there are larger voids and 
particles can move more freely. Let us notice that as far as the collective diffusion is involved 
the spread of a group of particles is faster in the presence of repulsive interactions than 
without. 46 

VI. CONCLUSIONS 

In this paper we have studied a one-dimensional system of inelastic hard-rods coupled 
to a stochastic heat-bath with the idea that it can represent a reference system in the area 
of granular gases to test theories and approximations. Due to the relative simplicity of the 
one dimensional geometry we have shown that it is possible to obtain relations between 
the macroscopic control parameters such as kinetic temperature, pressure and density. We 
tested these analytical predictions against the numerical measurements and found a fairly 
good agreement. It also appears that many properties of the heated one-dimensional inelastic 
hard rod system are similar to those of ordinary fluids. However, when we have considered 
how various physical observables fluctuate about their equilibrium values, many relevant 
differences have emerged. These range from the non-Gaussian behavior of the velocity 
distribution, the peculiar form of the distribution of distances between particles and of the 
energy fluctuations to the shape of the structure factor at small wave-vectors. Finally, we 
have found that the diffusive properties of the system are affected by the inelasticity and in 
particular the self-diffusion is enhanced. 

To conclude, in spite of the similarity between ordinary fluids and granular fluids, which 
has been recognized for many years and has made possible to formulate hydrodynamical 
equations for granular media in rapid, dilute flow, the presence of anomalous fluctuations 
in the inelastic case indicates the necessity of a treatment which incorporates in a proper 
way both the local effects such as the excluded volume constraint and the long ranged 
velocity and density correlations. Such a program has been partially carried out by Ernst 
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and coworkers^!, but needs to be completed regarding the description of the fluid structure. 



VII. APPENDIX A 



In the case of over-damped dynamics (i.e. large values of 7) one finds that the collective 
diffusion is given by 4 ^ 

(34) 



dp(x,t) = ld_( dp(p(x)) 
dt 7 dx \ dx 



where p(p(x)) is the local chemical potential. Expanding /i about its average value po we 
obtain: 

dp(p(x)) _ \5p] dp{x) ^ 



dx 

Substituting into Eq. (J3"4"j) we find 

dp(x,t) 1 



Sp 



dp(x) 



and using 



at 



S(0) 



7 



-Po 



5p 

Sp 
8p. 



dx 



d 2 p(x) 
dx 2 



K B T b 
Po 



dpo 



in the case of elastic hard rods we obtain: 



(36) 



(37) 



dp(x,t) lK B T b d 2 p(x) D d 2 p(x) 



dt 7 S(0) dx 2 5(0) dx 2 
Thus the renormalized diffusion coefficient is D = Dq/S(0). 



(38) 
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